Strong correlation induced charge localization in antiferromagnets 
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The fate of an injected hole in a Mott antiferromagnet is an outstanding issue of strongly correlated physics. 
It provides important insights into doped Mott insulators closely related to high-temperature superconductivity 
in cuprates |1-3|. Here, we report a systematic numerical study based on the density matrix renormalization 
group (DMRG). It reveals a remarkable novelty and surprise for the single hole's motion in otherwise well- 
understood Mott insulators. Specifically, we find that the charge of the hole is self-localized by a novel quantum 
interference mechanism purely of strong correlation origin |4, 5|, in contrast to Anderson localization due to 
disorders |6|. The common belief of quasiparticle picture is invalidated by the charge localization concomitant 
with spin-charge separation: the spin of the doped hole is found to remain a mobile object. Our findings unveil 
a new paradigm for doped Mott insulators that emerges already in the simplest single hole case. 



A critical step in understanding the physics of doped Mott 
insulators is to investigate the motion of a single hole in Mott 
antiferromagnets II2H51 171 1191 . It is easy to imagine that an 
injected hole in an antiferromagnetically ordered state leaves 
on its path the trace of spin mismatches in the staggered mag- 
netization direction (to be called "S'^-string" below). Such 
5^-string was expected |20| to cause the localization of the 
hole, but, as was later realized |9|, it can be self-healed via 
quantum spin flips. It was then conjectured [8l|9]|TTl[T2J that a 
quasiparticle picture should be still valid in long-distance, as if 
the Bloch theorem holds true even for such a many-body sys- 
tem. The quasiparticle picture was further supported by finite- 
size exact diagonalizations ||T3| on lattices up to 32 sites lfT4l . 
However, experimentally, in angle-resolved photoemission 
spectroscopy (ARPES) studies |21-23|, very broad single- 
particle spectral features have been observed in parent Mott 
insulator materials such as Ca2Cu02Cl2 and Sr2Cu02Cl2. 
They cast a serious doubt ifTSl [T8 1 on the rationale of un- 
derstanding an individual hole created by the photon within 
a quasiparticle description. Furthermore, the charge localiza- 
tion is a general phenomenon for a lightly doped cuprate be- 
fore the system becomes a high-Tc superconductor. Whether 
the localization here is intrinsic or due to the presence of dis- 
orders remains unsettled. 

Theoretically, the validity of a quasiparticle description has 
been questioned in the literature as well. It was first argued by 
Anderson ifTOl on general grounds that the quasiparticle spec- 
tral weight should vanish due to an unrenomalizable many- 
body phase shift induced when a hole is injected into a Mott 
insulator. Later, the discovery of the so-called phase string 
||4]|5l leads to a microscopic justification of this crucial obser- 
vation. Here, a phase string may be regarded as an 5* -string 
as opposed to the aforementioned "S"^ -string". It is irrepara- 
ble by quantum spin flips in contrast to the S'^-string. It was 
further predicted |18| that the phase string, instead of the S*^- 
string 1 20 1, is responsible for an intrinsic self-localization of 
the injected hole in two dimensions. 

In this work, a large-size numerical simulation is used as a 
powerful machinery to resolve this issue. This approach is of 
particular importance because the hole's motion is highly sen- 
sitive to quantum interference developing at long distance. To 
this end, we study the ladder systems by the DMRG method 
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Fig. 1: Quatum interference picture of one hole in the t-J model 
(A): As the hole moves from the injection point to a distant site, a 
sign sequence - the phase string |!4', '51 - is left behind as exempli- 
fied by the solid (green) and dashed (blue) lines. Here each + (— ) 
sign faithfully records an event that the hole exchanges its position 
with the nearest neighboring t (i) spin during its motion. (B) By 
contrast, the phase string can be wholly switched off in the o ■ t- 
J model (see the text) such that the amplitude associated with each 
path is positive definite (of the same absolute magnitude as in the t-J 
model). The presence (absence) of phase strings will drastically al- 
ter the long-wavelength behavior of the hole via quantum destructive 
(constructive) interference as demonstrated in this work. 



Il24ll . in which the length scale along one direction can be suf- 
ficiently large. As a ladder sample is long enough, we find that 
the charge is generally localized so long as the leg number is 
larger than one. The localization scale decreases monotoni- 
cally as the leg number increases, suggesting a stronger self- 
localization in the two dimensional limit. Contrary to this, if 
the sample is not long enough, the injected hole behaves itin- 
erantly, implying that the quasiparticle picture found in ear- 
lier numerical studies lfT3l [T4l is likely a small-size effect. 
We further show that complementary to the charge localiza- 
tion, the system exhibits novel spin-charge separation. That 
is, the "spinon" that carries the spin quantum of the doped 
hole remains mobile. We demonstrate that both charge self- 
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localization and spin-charge separation find their origin in the 
phase string, a unique property of the t-J model plfSllSSl. As 
a further evidence, we find that all of these exotic properties 
disappear once the phase string is turned off artificially in the 
simulation, and a well-defined quasiparticle description is re- 
covered in consistency with the common wisdom ||9l of spin 
polaron picture. 

Models and methods. — A prototypical doped Mott insula- 
tor is described by the t-J Hamiltonian, Ht-j ^ Ht + Hj, 
with 



(1) 



(iJ) 



Here, c]^. is the electron creation operator at site i, S; the 
spin operator, and ni the number operator. The summation is 
over all the nearest-neighbors, (ij) . The Hilbert space is con- 
strained by the no-double-occupancy condition, i.e., Ui < 1. 
At half-filling, Ui = 1, the system reduces to Mott insulators 
(antiferromagnets) with a superexchange coupling, J. Upon 
doping a hole into this system, J^i rii — N —1 (N the number 
of the lattice sites), and the hopping process is triggered as de- 
scribed by the hopping term, Ht, with t the hopping integral. 

For the single hole doped t-J model in a bipartite lattice, 
there exists an exact theorem, which states U |5] |25]| that the 
propagation of the hole is a superposition of quantum ampli- 
tudes of all the paths, with each path carrying a unique sign 
sequence known as the phase string, i.e.. 



(+1) X (-1) X (-1) X 



(2) 



Here, the sign ± on the right-hand side keeps track of an t 
or 4,-spin exchanged with the hole at each step of hopping as 
illustrated in Fig.[T](A). 

The quantum interference of phase strings from different 
paths is expected to be strong and play the fundamental role 
in dictating the long-distance behavior of the hole. In order to 
uniquely pin down the effect of phase strings, we also intro- 
duce the so-called a - t-J model with the hopping term Ht in 
Eq. ([T]i replaced by 



(3) 



(i.e., an extra sign a = ±1 is added). It is easy to show, fol- 
lowing Refs. im |5] |25]| . that the phase string (|2]i disappears, 
whereas the positive weight for each path remains unchanged 
[as illustrated in Fig. [T](B)]. In other words, one can artifi- 
cially switch on and off the phase string effect between the 
t-J and a ■ t-J models to study its novel consequences. 

Previously the DMRG algorithm has akeady been used to 
study the t-J model in ladder systems at low hole doping 
ifTTl l26l. Below we shall focus on the one hole case on bi- 
partite lattices of = x Ny, where and Ny are the 
site numbers in the x and y directions, respectively. By using 
the DMRG method, we shall study the ladders with small Ny 



(from 1 to 5) and sufficiently large N^- Here we set J as the 
unit of energy and focus on the t/J — 3 case unless otherwise 
specifically stated. Enough numbers of states in each DMRG 
block are kept to achieve a good convergence with total trun- 
cation error of the order of 10^^ — 10^^^. 

Self-localization of the charge. — Typical examples of the 
hole density distribution, (n^) = l — are shown in Fig.|2] 
(A) and (B) for A^j, = 3 and 4, respectively. Here the charge 
profiles are plotted along the x direction for a middle leg of the 
ladders. They are localized in the central region of the ladders 
with open boundaries. Upon summing up the distribution at 
all the sites of different legs, the sum rule: ^ ■ (n'') = 1 is 
satisfied. Examples of the contour plot of (n,^) in the x-y 



plane can be found in Fig. S3 in Supplementary Information. 

One may systematically change the ladder length to 
study the full width at half-maximum (FWHM) of the charge 
profile. Fig. [2] (C) shows that for ladders with different 
Ny, the FWHM increases linearly for small sample lengths, 
while saturates for sufficiently large sample lengths (except 
the Ny = 1 case). The former indicates that the doped hole 
still remains itinerant at a smaller A^^,. The latter suggests that 
a self-localization takes place for the doped hole, when the 
ladder length is long enough where the FWHM is no longer 
sensitive to the boundaries along the x-direction. Later on, 
by a more precise analysis, it will be shown that such local- 
ization only involves the charge of the doped hole, not its 
spin. Furthermore, Fig. |2] (C) clearly shows that the satu- 
rated FWHM at A'j^ > 1 monotonically decreases with the 
increase of the leg number Ny, implying even stronger local- 
ization in the two-dimensional limit. By contrast, there is no 
indication of saturation in FWHM for long one-dimensional 
chains (A^j, — 1, with N^. up to 500), which is consistent with 
the fact that the doped hole in strictly one dimension exhibits 
the Luttinger liquid behavior Q ■ 

The destructive quantum interference effect of phase strings 
in the t-J model provides a natural explanation for the charge 
"self-locaUzation" observed in the above DMRG results. To 
further confirm the phase string origin of self-localization, we 
switch off the phase string interference by two methods and 
then repeat the previous procedure of numerical simulations. 
In the first method, we study the a ■ t-J model introduced 
above. In the second method, we set the interchain hopping 
coefficient, t±, to be zero such that the hole can only move in 
the a; -direction as if in a one-dimensional chain. In both cases, 
the interference picture shown in Fig. [T] (A) breaks down. 
The numerical results for these two models are presented in 
Fig.|2](D) for two- and three-leg cases, in which the FWHM 
of the hole density distribution restores the behavior of linear 
increase with the sample length. It indicates that the charge 
becomes delocalized in both cases, in sharp contrast to the 
self-localization in the t-J model. 

Although the observed self-localization is insensitive to the 
parity (even-odd) of the leg number, the hole distribution is, 
as shown in Fig.|2](A) and (B) and Supplementary Informa- 
tion: for the even-leg ladders (A'j, ~ 2, 4), there are always 
small spatial oscillations on top of the Gaussian density pro- 
file, while they are absent for the odd-leg ladders (Ny = 3, 5). 
This can be attributed to distinct decaying behavior of spin- 
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Fig. 2: Self-localization of the charge in a 3-leg ladder system of size A'^ = 200 x 3 (A) and 4-leg ladder oi N = 200 x 4 (B). The density 
distribution (n'^) along the x axis in a middle leg is well fitted by the Gaussian function with a full width at half-maximum (FWHM). (C) 
The behavior of FWHM depends on the leg number Ny: for the ladders of Ny > 1, the FWHM first increases linearly at small Nx and then 
saturates at large Nx, with a length scale monotonically decreasing with the leg number; for A^j, = 1, the FWHM grows unboundedly as the 
sample length increases. (D) Delocalization of the charge in the absence of phase string interference for the two modified models, namely the 
a ■ t-J model and the t-J model with interchain hopping coefficient t± = 0. The corresponding FWHM unboundedly increases with Nx 
(open triangles, open diamonds, and crosses). The dashed line (= Nx/2) is a guide to the eye. 



spin correlations for the odd- and even-leg ladders at half- 
filling: the former (latter) follows a power (an exponential) 
law reflecting the absence (presence) of spin gap, see Supple- 
mentary Information for further explanations. Furthermore, 
the distinction of the even-odd effect disappears with the in- 
crease of t/J ratio (see Fig. S4i, indicating that the detailed 



spin dynamics, governed by J, is not essential to the charge 
self-localization effect. 

Momentum distribution. — Now let us examine the momen- 
tum distribution of the hole by studying n(k) = {c\^^c\^cj) , 
which can be obtained by a Fourier transformation of 
^^(c|^Cj(j). At half-filling, one finds n(k) = 1 and for the 
one hole case, 1 — n(k) is shown in Fig.|3] 

The insets of Fig.|3](A) and (B) present the hole momentum 
distribution 1 — n(k) as a function of for fixed ky = 27r/3 
in the 3-leg ladder and ky = 7r/2 in the 4-leg ladder, respec- 
tively. The value of ky is chosen in a way that the "sudden 
change" in 1 — n(k) can reach maxima by varying k^, accord- 
ing to the contour plots in the k^-ky plane (see Supplementary 
Information). A very interesting feature is that after the rescal- 
ing: 1 — n(k) [1 — n{k)]N, all the curves in the inset of 



Fig.|3jA) [or (B)] collapse onto a universal curve shown in the 
corresponding main panel. If one defines the Fermi surface by 
the sudden jump in the momentum distribution function, then 
the two universal curves suggest that the quasiparticle spectral 
weight is upper-bounded by 1 /N for large N, and vanishes in 
the thermodynamic limit. 

To satisfy the sum rule, X]k(l ^ '^(k)) = 1, the width of 
the jump in 1 — n(k) must remain finite in the limit iV — > oo. 
This is clearly shown in Fig. [3] (A) and (B), consistent with 
a finite localization length in the real space. We have also 
calculated the hole momentum distribution at different ratios 
oft/J. As shown in Supplementary Information, for a given 
sample size, the jump near the Fermi point is continuously 
reduced with increasing t/ J, which is qualitatively consistent 
with earlier work [13 ]. 

Again a sharp contrast arises once the phase string is turned 
off in the a -t-J model, as shown in Fig. [3jC) for the 4-leg lad- 
der Here the wide spread of 1 — n(k) in the momentum space 
seen in (A) and (B) collapses into a sharp peak at kx,ky = 0. 
In particular, a universal curve (the inset) is obtained by the 



rescaling: kr^ 



insteadof by 1 — [l~n(k)]N 
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Fig. 3: Momentum distribution. For both 3-leg (A) and 4-leg (B) 
ladders, the momentum distribution of the hole exhibits scaling be- 
havior: after the rescaling, 1 — n(k) — >■ [1 — n{'k)]N, the curves at 
different sample lengths in the inset collapse onto a single one repre- 
sented in the corresponding main panel. Note that we fix ky = ^ 
for the 3-leg case and ky = ^ for the 4-leg case, and plot the distri- 
bution along kx. Representive contour plots in the whole k^-ky plane 
are shown in Fig. |S5| (Supplementary Information). (C) The quasi- 
particle picture is recovered for the a ■ t-J model. The main panel 
shows that the momentum distributions of different sample lengths 
collapse to a single momentum at k^ — ky = for the 4-leg ladder; 
The inset shows that a universal curve is obtained for 1 — ?^(k) after 
the rescaling: kx — >■ N^kx, which suggests a nonvanishing quasipar- 
ticle spectral weight peaked at k = 0. 



as done in the main panel. This clearly indicates that the 
quasiparticle spectral weight is finite and a Bloch-like state 
is restored with a definite momentum at k = (a peak at the 
equivalent of kx,ky — tt is also found). 

Spin-charge separation. — Complementary to the charge lo- 
calization illustrated in Fig. |2] we further explore the fate of 
the spin-1/2 associated with the hole injected into the half- 
filled spin-singlet background. 

For the 3-leg ladder, the existence of spin-charge separation 
is explicitly demonstrated by Fig.|4jA). Here the single hole is 
created by removing a I spin electron. While the charge of the 
doped hole is well localized at the sample center, the extra spin 
of S*^ — 1/2 spreads over the sample. It is indicated by spin 
density (Sf) as well as the coarse grained (Sf)^ ^ , which is 
approximately uniform across the ladder This is probably the 
most direct display of electron fractionalization for a doped 
Mott insulator, thanks to the charge localization. Note that it 
is very different from the well-known one-dimensional Lut- 
tinger liquid, where both charge (holon) and spin (spinon) are 
mobile 121. 

For an even-leg ladder, the presence of a spin gap in the 
spin background renders the spin-charge separation less ex- 
plicit because the extra spin of S*^ = 1 /2 will also tend to stay 
in the hole region in order to reduce its superexchange en- 
ergy cost, as illustrated in Fig.QB) for the 2-leg ladder case. 
To directly observe a similar spin-charge separation picture 
in an even-leg ladder, the spin gap has to be small enough, a 
condition that can be achieved only if the leg number is suffi- 
ciently large. In this limit, the behavior of the single hole in 
both even- and odd-leg ladders is expected to converge with a 
stronger charge localization. 

Nevertheless, the spin-charge separation in the 2-leg ladder 
can be well confirmed by a different method. Note that, first, if 
the charge part - namely the holon - is localized, the localiza- 
tion center can be anywhere in a translationally invariant sys- 
tem. In particular, an "extended" profile as a linear superposi- 
tion of the localized wave packets at different positions should 
be also energetically degenerate. Second, if Fig. |4|B) really 
corresponds to spin-charge separation, the extra spin-1/2 - 
namely the spinon - should be more or less "free" within the 
regime expanded by the charge distribution. Then the spinon 
will prefer a more extended distribution of the hole to further 
lower its energy. In other words, the degeneracy of localized 
states should be lifted in favor of an extended profile in order 
to gain additional spinon energy. 

Indeed, by increasing DMRG sweep number with the trun- 
cation error reaching much lower than 10^**, the localized pro- 
file in Fig.|4]^B) is found to slowly converge to an extended one 
as shown in (C) with the total energy lowered by ~' 10^^ J. To 
identify the nature of this small energy gain, one may vary the 
sample size from — 16 to 800, and repeat the same pro- 
cedure to obtain the one-hole energy at each N^^ [with 
the same full convergence to the extended charge profile as in 
Fig.|4|C)]. As shown in the inset of Fig.|5jA), with a constant 
term excluded, E^^"^'' is well fitted by a{7r/Nx)'^ (blue trian- 
gles) with the coefficient a — 0.87 J, as if it is contributed by 
some "free particle" in a box of length iV^; 1271 . 

Next, we show that such l/N^ (finite-size) contribution to 
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Fig. 4: (A) Spin-charge separation in a 3-leg ladder with a J, spin electron removed. The charge of the hole is localized at the sample center and 
the spin of = 1/2 spreads over the entire sample, with the coarse grained density {Si)c.g. approximately uniform. (B) In a 2-leg ladder, 
the charge and spin of the doped hole remain in the same localized region because of a finite spin-gap in the spin background. (C) An extended 
profile of the charge can be obtained with further reducing DMRG truncation error (see the text); In this case, although the charge remains 
localized, the mobile spinon gains a tiny additional energy (cf. Fig.|5]l. (D) In a = 400 x 2 ladder, by adding weak chemical potentials the 
hole density can be localized anywhere along the ladder. 



^i^-hoie jjQj. fj-Qjjj j-jjg charge sector. For this purpose, we 
make the 2-leg ladder as a ribbon by connecting the bound- 
aries on the two edges of the ladder in the x-direction. Then 
we calculate the energy difference AEq'°^'^ = Eq^°^''{'k) — 
^^hoie^Q-j ^jjjj threading through a flux tt into the ribbon. It 
corresponds to the change of boundary condition from the pe- 
riodic to anti-periodic one for the hopping term and is equiv- 
alent to coupling a small electric field to the charge sector. As 
shown in Fig. D[ A), AEq^°^'^ oscillates strongly and falls off 
exponentially as e^^^^^ with ^ ~ 14.5. It behaves completely 
different from 1 /N^ behavior exhibited in Eq^°^'', indicating 
that the charge indeed remains localized. The "free particle- 
like" behavior in should therefore solely originate from 
the charge-neutral part of the doped hole, i.e., the spinon. 

By a sharp contrast, upon turning off the phase-string in the 
a-t-J model, both £'^-h°''= and AS^*°''= follow the same l/iVj 
behavior as clearly shown in the inset of Fig.|5]^A). It is con- 
sistent with the previous conclusion that in this case the doped 
hole is restored as a quasiparticle, which carries both charge 
and spin and thus responds very similarly to the boundary con- 
dition (i.e., the finite-size effect) as a single mobile object. 

The spin-charge separation is well-known for the one- 
dimensional chain |2|. Figure |5jB) presents the correspond- 
ing 1-leg results. As shown in the inset, the finite-size effect of 
i?^*"'*^ (cx is indistinguishable from the phase string 



free case (with open boundary conditions). But the charge re- 
sponse of AEq^°^'' is again totally different [the main panel of 
Fig.[5jB)]. It oscillates strongly as a function of A^^^, similar 
to the 2-Ieg case in Fig.|5jA). As a matter of fact, it can be di- 
rectly attributed to the phase string modulation in the expres- 
sion for the ground state energy given in Ref. 1 18]. The main 
difference from the 2-Ieg case is that the envelop of AEq^°^'' 
decays with increasing in a power law fashion (cx 1/-/VJ) 
instead of an exponential one. This is consistent with the fact 
that the holon is not localized in the one-dimensional case. 
It is interesting to notice that even without the phase string 
effect, £'^-h°''= and AE:',*'"'" behave distinctly for the a ■ t-J 
model. The former cx l/Nl'^, while the latter cx 1/A^J, also 
indicating a spin-charge separation although the distinction is 
much less dramatic as compared to the t-J model. 

The presence of spin-charge separation enriches substan- 
tially the behavior of the hole localized states in the t-J ladder 
systems. Due to the residual spin-charge interaction, the de- 
generate charge localized states are superposed to form a more 
"extended " distribution to minimize the spinon energy. It has 
been shown |28 | that generally the DMRG method tends to 
pick up a minimally entangled state among quasi-degenerate 
states. Thus, spatially localized ones such as those in Fig. [2] 
are naturally obtained because a further spinon energy gain 
for a bigger FWHM would be too tiny to be detected with 
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Fig. 5: The energy difference AEq^"^'' between anti-periodic and pe- 
riodic boundary conditions along the x-direction for the charge shows 
an oscillation vs. Nx, which is a direct manifestation of the phase 
string effect. The exponential decay of the envelop (dashed curves) 
of AEq^^"'^ shows the charge localization for the 2-leg ladder (A), 
while the power law decay for the 1-leg ladder case (B) indicates a 
Luttinger liquid behavior. By sharp contrast, in the inset of (A), a 
completely different finite-size behavior (l/N^ ) is exhibited in the 
ground state energy Eq""'" for the 2-leg ladder, confirming the spin- 
charge separation. For comparison, both AEq^"^'' and Eq^"^'' follow 
the same l/N^ dependence once the phase string is turned off, as 
shown in the inset of (A), where the quasiparticle picture is restored 
with a spin-charge recombination. In the inset of (B), similar com- 
parisons for the 1-leg case are presented. 



a given small DMRG truncation error. It also explains why 
the FWHM of charge distribution shown in Figs. |2]and|4]is 
usually bigger than the intrinsic localization length £^ as ex- 
hibited by Ai?Q*°''^. But the former is still valid to charac- 
terize the general trend of the localization, say, as a function 
of the leg number in Fig. [2] On the other hand, by adding 
some weak local chemical potentials, the quasi-degeneracy of 
the hole states can be lifted such that the hole density profile 



is truly localized around the impurities at different locations 
[as indicated by different colored curves in Fig. |4|D) for a 
N = 400 X 2 ladder]. Here the increase of the total energy 
of the original t-J Hamiltonian is weak (about 10^* J), pre- 
sumably from the spinon part as the FWHM is still larger than 

e 

Discussions. — Phase string, charge localization, and spin- 
charge separation are striking results of a single-hole-doped 
Mott insulator. They are truly incompatible with a quasipar- 
ticle description. However, for their novel effects to get fully 
unveiled, a large sample is needed. Instead of tackling directly 
a two-dimensional square lattice, in this work, we have chosen 
ladders that DMRG can access. As it turns out, very surpris- 
ing consequences do show up so long as the ladders are long 
enough, and, there are more than one paths to realize quantum 
interference, e.g., in a two-leg system. 

The most essential physics underlying the observed exotic 
phenomena is that the motion of the single hole acquires a se- 
quence of signs that precisely and permanently record each 
microscopic hole-hopping-spin-backflow event in the quan- 
tum spin background. These irreparable phase strings cause 
singular and destructive interference of quantum waves of the 
charge. The longer the paths are, the stronger the destruc- 
tive interference is. The present results show unequivocally 
that the self-localization of the charge is a combined effect 
of quantum interference and multiple scattering between the 
doped charge and spin background, in analogy to Anderson 
localization |6|. Contrary to this, the charge remains delo- 
calized in the single chain case due to the lack of quantum 
interference involving different paths. 

The charge localization is further found to be concomitant 
with spin-charge separation for the two and more-leg ladders. 
Whereas the phase string renders kinetic energy suppression 
of the charge, the spin associated with the doped hole remains 
unfrustrated and thus maintains a coherent motion as a spinon. 
This has been clearly demonstrated by the distinct finite-size 
effects of the charge and spin degrees of freedom. In this 
sense, the electron fractionalization found in the single hole 
case can be attributed to the energetic reason. Note that the 
cause for charge localization, i.e., the nontrivial phase string 
effect, requires an uncertainty of the spin quantum associated 
with the spin backflow (leading to the ± signs). If each spin 
backflow could maintain a definite spin polarization direction, 
the phase string effect would become trivial. Simultaneously 
the hole would also carry a definite spin as a quasiparticle. 
But such a spin-charge recombination, being extremely tight 
in space, would cost too much kinetic energy. 

This is fundamentally different from the premise for a spin 
polaron picture - a conventional wisdom in which the distur- 
bance of the hole hopping to the spin background is described 
by virtual spin excitations that eventually relax back after the 
hole moves away. In fact, the spin-polaron or quasiparticle be- 
havior has been shown to be recovered once the phase string 
is "switched off in the so-called a ■ t-J model, where every- 
thing seems "normal" in a conventional sense. Here a spin- 
charge recombination (not too tight spatially) with a good ki- 
netic energy for the quasiparticle can be realized since there is 
no more singular phase string effect to strongly suppress the 
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kinetic energy of the charge. 

Our results have far-reaching consequences. First of all, 
they suggest that a new state of quantum matter may natu- 
rally emerge in sufficiently low doping. Such a state is insu- 
lating because of the self-localization of doped charges, not 
because of a charge gap like in the undoped case. This is a 
novel non-Anderson-type localization phenomenon since ex- 
ternal disorders are absent. Secondly, the strong suppression 
of the kinetic energy results in a large quasi-degeneracy in 
the low-energy states. It is expected to lead to a "glassy" 
behavior in the presence of weak external impurities and a 
"competing-order" phenomenon in the presence of additional 
weak interactions. It also provides a possible explanation for 
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Supplementary Information 



In this supplemental material, we present some technical details. 

Spin-spin correlations at half-filling. — At half-filling, the t-J model reduces to the Heisenberg model. For the isotropic 
Heisenberg coupled-chain systems, the behavior of the even-leg ladders is dramatically different from that of the odd ones. A 
well-known fact is that the even-leg ladders have a spin gap while the odd-leg ladders are gapless. The former leads to an 
exponential decay of the spin-spin correlation while the latter to a power law decay (see Fig.jSTJi. The spin gap for the even-leg 
ladders is expected to vanish in the large leg number (namely two-dimensional) limit. Indeed, the spin structure factor in the 
4-leg case has already exhibited strong antiferromagnetic correlations as shown in Fig. S2 These results are consistent with 
earlier DMRG work|30 |. 

Some details of the hole density distribution. — Fig. S3 presents the contour plot of (n^) for different ladders, iV — 40 x Ny 
with Ny = 2, 3, 4, and 5. For each contour plot, we have checked numerically that the sum rule is satisfied, i.e., '}2n {n^) = 1, 
where the summation is over all the sites. 



A prominent feature can be clearly seen in these contour plots (see also Fig. S4 (B)). That is, for even-leg ladders there is 
a spatial oscillation on top of the localized density profile, while it is absent for odd-leg ladders. We wish to point out that 
this is concomitant with the parity effect of the spin-spin correlation shown in Fig. SI In fact, the spin-gap effect is reduced 
as the hopping integral t increases, and eventually the parity effect disappears in the large t/ J limit. Accordingly, the spatial 
oscillation of (nj*) must be diminished also in this limit. This prediction has been fully confirmed by numerical simulations, see 



Fig. S4 (B). In addition, we have observed that the localization length monotonically decreases as the ratio t/ J increases, see 
the insets of Fig. S4 This suggests that the spin dynamics, governed by the superexchange J, does not play important roles in 
establishing the observed self-localization. 

Electron and hole momentum distributions. — Fig. S5 shows the contour plot of the hole momentum distribution function 
1 — n(k) for the 3-leg (A) and 4-leg (B) ladders, respectively at — 80. The maximum of 1 — n(k) appears at ky=zLTr/2 
for the 4-leg ladder and fcy=±27r/3 for the 3-leg ladder The </ J-dependence of the hole momentum distribution is shown in 



Fig. S6 They show that the momentum distribution jump is monotonically reduced with the increase oft/J. 



[30] S.R.White, R. M. Noack, and D. J. Scalapino, Resonating Valence Bond Theory of Coupled Heisenberg Chains. Phys. Rev. Lett. 73, 
886(1994). 




5 10 15 20 25 30 5 10 15 20 25 
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Fig. S 1 : The main panels are the spin-spin correlation | (^S'f 5| ) | in the 3-leg (A) and 4-leg (B) ladder, respectively, ix and jx are chosen to be 
on the middle leg of the corresponding ladder. The insets are In-ln plot (A) and semi-ln plot (B). They are well fitted by straight (blue) lines, 
indicating a power law decay for the 3-leg ladder and an exponential decay for the 4-leg ladder. 
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Fig. S2: Spin structure factor for the 4-leg ladder (A'^=80x 4). Tfie peak is located at (tt.tt), indicating strong spin antiferromagnetic correlations 
even though a finite spin gap is present. 




Fig. S3: Contour plots of the hole distribution in real space. The sample size 40 x 2 (A), 40 x 3 (B), 40 x 4 (C), and 40 x 5 (D) from top 
down. 
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Fig. S4: The hole density profile (nj") (with i located on the middle chain) at different ratios of t/J. The ladder size is 80 x 3 (A) and 80 x 4 
(B). The localization length (the FWHM) is found to decrease as t/J increases (insets). Moreover, the amplitude of the spatial oscillation in 
the 4-leg ladder is suppressed as t/J increases, indicating that the behavior of even- and odd-leg ladders converges in the large t/J limit. 
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Fig. S5: Contour plots of the hole momentum distribution 1 — n(k) for 3-leg (A) and 4-leg (B) ladders. 




Fig. S6: The hole momentum distribution at different ratios oft/ J for 3-Ieg (A) and 4-leg (B) ladders. 



